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Abstract 

We review a relativistic approach to the heavy quark physics in lattice QCD by 
applying a relativistic 0(a) improvement to the massive Wilson quark action on the 
lattice. After explaining how power corrections of niqa can be avoided and remaining 
uncertainties are reduced to be of order (otIqcd)^, we demonstrate a determination of 
four improvement coefficients in the action up to one-loop level in a mass dependent 
way. We also show a perturbative determination of mass dependent renormalization 
factors and 0(a) improvement coefficients for the vector and axial vector currents. 
Some preliminary results of numerical simulations are also presented. 
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§1. Introduction 



Lattice QCD should allow quantitative predictions for the heavy quark physics from first 

principles. Up to now, however, most approaches have based on the nonrelativistic effective 

theory, with which the continuum limit can not be takenP In this report we review a new 

relativistic approach to the heavy quark physics in lattice QCD developed in a series of 
publicationsPEJaJS)E)EHH) 

We consider a relativistic 0(a) improvement to deal with the heavy quarks on the lat- 
ticed We discuss cutoff effects by extending the on-shell improvement programm^^''^'^''^ 
from massless to massive case. An important finding is that in our formulation leading 
cutoff effects of order {rriQa)" are absorbed in the definition of renormalization factors for 
the quark mass and the wave function. After removing the next-leading cutoff effects of 
0((mQa)"a/lQCD) with four parameters in the quark action properly adjusted in a mgo de- 
pendent way, we are left with at most ©((o/Iqcd)^) errors. We show a determination of 
the four parameters in the quark action up to one-loop level for various improved gauge ac- 
tions.!^ They are determined free from the infrared divergences, once their tree level values 
are correctly tuned in the mga dependent way. 

We also make the 0{a) improvement of the vector and axial vector currents at the 
one-loop level.!^ We give a general discussion about what kind of improvement operators 
are required from the symmetries allowed on the lattice, in which the Euclidean rotational 
symmetry is violated because of mga corrections. We consider both the heavy-heavy and 
heavy-light cases. 

This report is organized as follows. In Sec. |21 we consider a relativistic 0{a) improve- 
ment to handle the heavy quarks on the lattice avoiding large mqa corrections. The four 
parameters in the quark action are determined up to one-loop level in Sec. 0] In Sec. E]we 
explicitly show the 0(a) improvement of the axial vector current up to one-loop level. With 
the use of the 0(a) improved quark action and axial vector current, we make some numerical 
studies in quenched QCD focusing on restoration of the space-time symmetry. Their results 
are presented in Sec. El Our conclusions are summarized in Sec. EJ 

§2. On-shell improvement of the massive Wilson quark action 

We consider a relativistic 0(a) improvement to control mga corrections for the heavy 
quarks on the lattice. The basic idea is to apply the on-shell improvement programme 
which has been developed in the small mass case, to the heavy quarks on th lattice. This 
method allows us to obtain the physical quantities in the continuum limit without requiring 
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harsh condition niQa <^ 1 that is not achievable in near future. 

Before going into details, we first remark on the on-shell improvement. As explicitly 
stated in Ref . IT^ . the on-shell improvement is meant to improve the correlation functions in 
which local composite fields are separated by non-zero physical distances. All on-shell quanti- 
ties (particle energies, scattering amplitudes, normalized matrix elements of local composite 
fields between particle states, etc.) are extracted from these correlation functions. It 
should be stressed that the on-shell quantities are not restricted to the spectral ones. 

Let us consider general cutoff effects for the heavy quarks on the lattice, where the 
heavy quark mass rriQ is allowed to be much heavier than Aqqd- Under the condition that 
niQ ^ Aqcb and mga ~ 0(1), we assume that the leading cutoff effects are 



where fi{mQa) {i > 0) are smooth and continuous all over the range of mqa and have Taylor 
expansions at mqa = with sufficiently large convergence radii beyond mqa = 1, taking 
/o(0) = and /i(0) ~ 0(1) for i > 1. The essential point in this assumption on cutoff effects 
is that the cutoff effects of 0(a/lQCD) in the chiral limit are still 0{aAQCB) even if the quark 
mass is increased. This means that our power counting is not based on the nonrelativistic 
effective theory. To control the scaling violation effects we want to remove the cutoff effects 
up to fi{mQa)aAQCD by adding the counter terms to the lattice quark action with the on-shell 
improvement. If mqa is small enough, the remaining /2(mQa)(aylQCD)^ contributions can 
be removed by extrapolating the numerical data at several lattice spacings to the continuum 
limit. Otherwise, in case of sufficiently small lattice spacing, the 0((a/lQCD)^) errors can be 
neglected. 

We first search for the relevant counter terms required in the on-shell improvement. Listed 
below are the allowed operators under the requirement of the gauge, axis interchange and 
other various discrete symmetries on the lattice, where the chiral symmetry is not imposed. 
According to the work of Ref. EJ, all the operators with dimension up to six are given by 



foimqa) > fi{mQa)aAQCD > f2{mQa){aAQCDf > ■ ■ , 



(24) 




(2-2 






(2-4 






(2-6; 

(2-7; 
(2-8; 
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06d(x) = iq{x)-if,[D„ F^,]q{x), (2-9) 
Oeeix) = q{x)p\{x), (2-10) 
0^f{x) = q{x)rq{x)q{x)rq{x), (2-11) 

where -T = 1, 75, 7^^, 7^75, cr^i/- These operators lead to a following generic form of the quark 
action on the isotropic lattice: 

•^f' = 5Z C3O3(a:) + C404(a;) + 5^C5i(95.(x)+ ^ cM,{x) + ■ ■ ■ , (2-12) 

X L i=a,b i=a,...,f 

where C3, . . . ,C6/ are functions of the bare gauge coupling g and the power corrections of 
rriQa. At mq = the contributions of dimension 6 operators are of order (o/Iqcd)^, which 
are negligible for the 0{a) improvement. However, we are now interested in the case of 
mq > Aqcd and mga ~ 0(1). 

We first point out that, regardless of the magnitude of mqa, the mqa corrections to the 
quark mass term and the kinetic term can be absorbed in the renormalizations of the quark 
mass and the wave function Zg. For the sake of convenience we choose C3 = mo and 
C4 = 1. 

In the next step we reduce the number of basis operators with the aid of the classical 
field equations. It is easily found that O^a, C'efe, C'ec and Oge can be related to the quark 
mass term or the kinetic term. In the on-shell improvement these operators are redundant 
and can be eliminated from the action of eq. ()2-12|) . The operator O^a, however, is used to 
avoid the species doubling and the value of its coefficient c^a is given by hand. 

The remaining operators are O^h, C'ea, C^ecz and Og/- It is easy to see that Oq^ and 
O^f are ©((/iQCDCt)^); and therefore are irrelevant for the 0{a) improvement. The operator 
O^b is the so-called clover term, for which the nonperturbative method to determine the 
coefficient c^b in the massless limit is already established.!^ However, the contributions of 
[mqaYO^h in > 1) cannot be neglected in the present condition that allows mqa ~ 0(1). 
For 0(a/lQCD) improvement the coefficient csb has to be adjusted in the mass dependent 
way. 

Our main concern is Cga- Under the condition that 9og(x) ~ mqq{x) and diq{x) ~ Piq{x) 
with mq ^ \pi\ and mqa ~ 0(1), we have to treat the time and space components of O^a 
in a different way: They follow different power counting. Actually, the contribution of the 
time component q{x)'yoDQq{x) , which is 0((ylQCDa)^) in the massless limit, is not negligible 
any more in the present condition. It can be related to other lower dimensional operators 
by the equation of motion: 

a'^q{x)-foDlq{x) = --{mqa)^q{x)q{x) 
a 
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-{mQayq{x)-fiDiq{x) 
+a{mQa)q{x)D^q{x) + 0{{AQCBaf). 



(2-13) 



This relation tells us that the contribution of q{x)^oDQq{x) is expressed by the lower di- 
mensional operators multiplied by the power corrections of mqa, so that the coefficients of 
the space derivative terms in O4 and O^a become different from those of the time deriva- 
tive terms. On the other hand, the contribution of the space component of Oqq is to be 
©((ylqcDo)^), which is neghgible for the 0{a) improvement. It is essential to note that the 
Lorentz non-covariant terms like Oga yield the difference of magnitude between the time and 
space components due to finite rriQa corrections. 

The generalization of the above argument to any operators with higher dimensions makes 
the discussion more transparent. Let us consider an arbitrary operator with A + k dimension, 
a'^Oi+k, where we write the lattice spacing a explicitly. The operator O^+k contains I pairs 
of q and q and n covariant derivatives with A + k — 3xl + n. Using the classical field 
equation, some (but not all) of covariant derivatives can be replaced by the quark mass 
niQ. For I > 2 the largest possible power of the scaling violation is (mQa)"(aylQCD)^'~^- 
Therefore the operators which contain four or more quarks are irrelevant for the O(aylQCD) 
improvement. All the relevant contributions come from the quark bilinear operators. With 
the aid of the classical field equations, they can be reduced to 



for n > 0. The time and space components of O4 and 05a,5b should be treated separately 
in case of finite m^a, where the space-time asymmetry refiects the contributions of the 
higher dimensional operators that break the rotational symmetry. Now we know that the 
seven operators are needed for the O(aylQCD) improvement. Since three coefficients among 
these seven operators can be absorbed in Z^, Zg and the Wilson parameter for the time 
derivative of O^a as already explained, the remaining four coefficients have to be actually 
tuned. 

In conclusion, at all order of m^a, the generic quark action is written as 




(2-14) 
(2-15) 





(2-16) 




(2-17) 



= X] moq{x)q{x) + q{x)'yoDoq{x) + q{x)'yiDiq{x) 



X 
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2 .V / u.v . 2 . 

i 

ina 1 — 

2-11 



-^-^CE^q{x)aoiFoiq{x) - ^-^CB^q{x)aijFijq{x 



2 , V. .... , ^ 

where we are allowed to choose = 1 and the four parameters u, r^, ce and cb are to be 
adjusted. In general these parameters have the form that X = X*^"^((7^)(mQa)" with 
X = u, Ts, Ce and cb, and should agree with the one in the massless 0(a) improved 
theory: z/(°) = 1, ^ = n = 1, = = Csw^ Note that u = I + 0((mQa)2) and 
Ts = Tt + 0{mQa) since the space-time asymmetry arises from Lorentz non-covariant terms 
such as O^a via the on-shell reduction of eq. ()2-13p accompanied by power corrections of mqa. 
In brief the differences between u and 1, and r^, ce and cb reflect the contributions of 
Lorentz non-covariant terms with higher dimensions. 

We find that a similar quark action for the heavy quarks has been proposed in Ref. |[T5|) . 
The important difference is that the parameter is redundant in their formulation, while it 
should be tuned in ours. In this sense our action is equal to theirs with a special choice of 
parameters. The reason of discrepancy in the number of relevant operators between in our 
formulation and that in Ref. IT^ is explained in detail in Ref. We explicitly show in the 
next section that the parameter actually needs to be adjusted to reproduce the correct 
on-shell quark-quark scattering amplitude. 

§3. Determination of the improvement parameters in the quark action 

3.1. Tree level 

The four improvement parameters in the quark action of eq. ()2-18|l are determined such 
that 0(a) cutoff effects in on-shell quantities (particle energies, scattering amplitudes, nor- 
malized matrix elements of local composite fields between particle states, etc.) are removed. 
We employ the on-shell quark-quark scattering amplitude to determine the parameters u, rs, 
Ce and cb at tree level, which are adjusted to reproduce the continuum form of the scattering 
amplitude removing the mqa corrections,!^ 

T = -g'^{T^fu{p)-f^u{p)D^y{p - p')u{q')-i^u{q) 
-9^{T^fu{q)-i^u{p)D^y{p - q')u{p')-fM(l) 

+0{{p,a)Mq,ay,{p[a)Mq',af), (3-1) 

where p, q denote the incoming quark momenta and p', q' for the outgoing quark ones. D^^, 
denotes the gluon propagator. This improvement procedure follows the previous woriP' 
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that determined the csw = cb = ce parameter up to one-loop level in the massless case. At 
the tree level the quark-quark-gluon vertex is written as 



u{p')A^o\p,P>{p)) = (^(p')^7ou(p))_t + 0(fea)^ (p^a)^), (3-2) 

/ latt 

= (^(p')^7.n(p))cont + 0{{v.a)\ {p^af), (3-3) 

latt 



for 



.(0)/ /X . fPo+Po\ . . fPo+Po\ 

+ ^cos J]ao.sin(p. -pD, (3-4) 



^(0)/ >^ ■ (0) (Pk+P'k\ , (0) • fPk+P'k 

[P^P) = Ik COS I ^ 1 + r^^' sm I ^ 



fPk-p'k\ . I IS 

COS I — ^ — 1 afco sm(po - Pq) 

(0) / _ ' \ 

V\ cos ( ) sin(p, - pj), (3-5) 



where the spinor on the lattice is given by 



with A^(p) = (— ■i)sin(po) + + — cos(po)) + ~ The 0(a) improvement 

condition yield^ 

.(°) = (3.7) 

_ cosh(m^°^) + sinh(m^°^) sinh(mf'^) 
m), my 

cf=nv^'\ (3-9) 

4°)=rf, (3-10) 

where mp°^ is the tree-level pole mass explained below. It is now shown that the four 
parameters are uniquely determined from the on-shell quark-quark scattering amplitude. 
This is an evidence we actually need four improvement parameters in the quark action of 
eq.dHHl)- 
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It is instructive to show that the v and parameters are also determined from the quark 
propagator which is obtained by inverting the Wilson-Dirac operator in eq. ()2-18|) . 

^q^{v) = ^7osin(po) + 7ism{pi) + itlq 

i 

+rt{l - cos(po)) + ^(1 - cos(pi)), (3-11) 

i 

At the tree level the parameters are adjusted such that the above quark propagator repro- 
duces the correct relativistic formP 

SqVP) = —^^ 7^2 + (no pole terms) + 0((pia) ) (3- 12) 

around the pole. Z^q^ and m'p^ are extracted with = 0, 



mf^ = log 



'mo + rt + \/ml + 2rtmo 



) = cosh(mJ°)) + rtsinh(mj,°)). (3-14) 

Imposing finite spatial momenta the parameter u is determined by demanding the correct 
relativistic spinor structure on S^^^p) of eq. ()3-ll|) . Comparing the coefficients of 70 and 7^ 
in the numerator we obtain 

.(°) = ""^^5'^ - (3-15) 
rrip 

The parameter is determined such that the correct dispersion relation is reproduced: 

E' = ml + J2pl + 0ipt)- (3-16) 



The result is 



(0) _ cosh(m), ') + rtsinh(m), ') sinh(m), ') 

(Zf-z/W). (3-18) 



(0) 

mi, 



It should be noted that the values of z/(°^ and ri°^ are exactly the same as those determined 
from the on-shell quark-quark scattering amplitude. This is not an accident: The correct 
relativistic form of quark propagator yields the correct relativistic form of Dirac spinor 
required in the calculation of matrix elements. Since it is simpler to treat the conditions on 
the quark propagator, we employ them to determine the one-loop contributions to v and 
in the next subsection. 



8 



3.2. One-loop level 

The one-loop contributions to the quark self-energy come from the rainbow and the 
tadpole diagrams, which are written as 



i-fo smpoBo{p, mo) + z/i ^ 7^ sin piBi{p, mo) + C{p, mo) 



Incorporating this contribution, the inverse quark propagator up to the one-loop level is 
written as 

S~^{p, m) = i'jo sinpo[l - g'^Bo{p, m)] + z/i ^ 7^ sinpi[l - g'^Bi{p, m)] + m 

i 

+2nsm' (^1) +2r,5^sin2 (|) _/(7(p,m), (3-19) 

i 

where we redefine the quark mass as 

m = mo- g^C{p = 0,m = 0), (3-20) 
C{p,m) =C{p,m) -C{p = 0,m = 0). (3-21) 

With this definition the inverse quark propagator satisfies the on-shell condition for the 
massless quark up to the one-loop level : S~^{pq = 0,pj = 0, m = 0) = 0. 

The u and parameters are determined by employing the same improvement condition 
as the tree level. The parameter u is determined from the relativistic spinor structure in 
of eq. ()3-19|) at the pole. 

u[l-g'B,ip\m)] = '^^^^[l-g'Bo{p\m)], (3-22) 

nip 

where p* = (po = i'rrip^Pi = 0). We show the quark mass dependences of u^^y^^^^ over the 
range < m^^ < 10 for the plaquette and the Iwasaki gauge actionJ^' in Fig. ^a). The 
solid lines depict the results of the interpolation with a rational expression. The errors are 
within symbols. The parameter is determined from the relativistic dispersion relation as 
done at the tree level. The m^^ dependences of ri^^r^p for the plaquette and the Iwasaki 
gauge actions are shown in Fig. E^b). 

Let us turn to the one-loop calculation of ce and cb- Recently the authors have shown 
the validity of the conventional perturbative method, which uses the fictitious gluon mass 
as an infrared regulator, to determine the clover coefficient csw up to the one-loop level 
in the massless case from the on-shell quark-quark scattering amplitude.!^ We extend this 
calculation to the massive case. According to Ref . ITE|l . it is sufficient to improve each on-shell 
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Fig. 1. (a)z^(iVi^^°^ and {h)ri'^^ /A^^ function of mp*^ . Circles denote the plaquette gauge action 
and triangles for the Iwasaki gauge action. 



quark-quark-gluon vertex individually. To determine the one-loop coefficients and 
we need six types of diagrams shown in Fig. |21 We ffist consider to calculate c^j^\ Without 
the space-time symmetry the general form of the off-shell vertex function at the one-loop 
level is written as 

{P, q,m) = [p, q, m) 

i=a,...J 

+MF^ + MoF^ + iolk^F^ + -ik'^o'^F^ + (i(io7kF^o 

+{pk + Qk) [^f + + 4H^^ + ml] 

HPk - qk) [G\+iGl + iGl + iiGf\ + 0(a2), (3-23) 

where Ak{p,q,m) = /if ^(p, g, m)+5(M^^^(j9, g, m)+0(/) and^ = ELo^^^a, <i = ELo^^^a, 
^0 = Polo, io = qolo- The coefficients F^^ {i = 1,...,10), {i = 1,2,3,4) and 
{i = 1,2,3,4) are functions of p^, q"^, p ■ q and m. From the charge conjugation symmetry 
they have to satisfy the following condition: = F^, F^ = F^, F^ = F^, Fg = Ff^Q, 
H2 = H^, Gi = G'l = and G2 = — Cg. Sandwiching yl[.^'*(p, g, m) by the on-shell quark 
states u{p) and u{q), which satisfy ^m(p) = impu{p) and u{q)4 = impu{q), the matrix element 
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(a) 



(b) 






Fig. 2. Quark-quark-gluon vertices at one-loop level, p (q) is incoming (outgoing) quark momen- 
tum. 

is reduced to 



u{q)A'j}\p,q,m)u{p) 
u{q)lku{p) {Ff + imp{F!^ + F, 



+ {pk + qk)u{q)u{p) [H^ + imp{H!^ + H, 
+ {Pk - qk)u{q)u{p) {G\ + imp{G 



m. 



2TTk\ 

G^}+0(c 



(3-24) 



where we use 



Ft F!^ 



F| and F| 



(0) 



(Note that we can replace nip with 
in the one-loop diagrams.) The first term in the right hand side contributes to the 



renormalization factor of the quark-quark-gluon vertex, which is equal to z'^^ at the tree 
level. With the use oi G\ = G\ = and G^ = —G\, we find that the last term of eq. ()3-24|) 
vanishes: this term is not allowed from the charge conjugation symmetry. It is also possible 
to numerically check G\ + imp{G2 + G^) — rripG'l = 0. 

The relevant term for the determination of cb is the third one, which can be extracted 
by setting p = p* = (po = irrip, pi = 0) and q = q* = {qo = imp, qi = 0) in eq. ()3-23|) : 



in. 



i^k 



(3-25) 
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where we have used the fact that F^, and are functions of p^, and p ■ q, so that 



dpi 



dpi 



p=p',q=q-' 



p=p* ,q=q* 



dqi 



dpi 



dqi 



dqi 



p=p',q=q' 



p=p* ,q=q* 



p=p'',q=q'' 







(3-26) 

(3-27) 
(3-28) 



P=P ,9=9 

with j = 1, . . . , 10, / = 1, 2, 3, 4 and z = 1, 2, 3. 

We should remark that the third term in eq. ()3-24|) contains both the lattice artifact of 
0{pka, qko) and the physical contribution of 0{pk/m, qk/m). The parameter cb is determined 
to eliminate the lattice artifacts of 0{pka,qka): 



(1) 

- r, 



(1) 



I latt 

'p=p*,q=q* 

-Zf^ [Hi + rm,{H!^ + Hi) - m'^Ht] 



, . (3-29) 

where we take account of the tree-level expression for the quark-quark-gluon vertex in eq. ()3-3|) 
and eq.(3.51) of Ref. We show the quark mass dependences of c^^ I ^^b ^^^"^ range 
< mp°^ < 10 for the plaquette and the Iwasaki gauge actions in Fig. IHl^a). The solid lines 
represent the interpolation with a rational expression. The errors are within symbols. 

The calculation for c^^'' is done in a similar way as for c^^'': we determine c^^'* to remove 
the 0{a) contribution from the on-shell matrix element u{q)A\^\p,q,m)u{p). The quark 
mass dependences of c'e' / c'e' shown in Fig. IHfb). 

Before closing this section we should remark an important feature in the calculation of 



,(1) 



c]j' and c"q\ The infrared divergences originating from Figs. 121(a), (b), (c) contain both the 



lattice artifacts and the physical contributions. The former exactly cancels out if and only if 
the four parameters rf\ c^^ and are properly tuned. In Ref. E} we demonstrate it 
by explicitly writing down the infrared behaviors of the one-loop diagrams. This is another 
evidence that the tree- level improvement is correctly implemented in Ref. E|- It may be also 
instructive to mention the massless case,'^''^^ where we find a similar situation: The infrared 
divergences originating from Figs. |21 (a), (b), (c), (e), (f) cancel out if and only if csw*-"-* is 
correctly adjusted at the tree level, namely csw''"'' = 1- 

§4. 0{d) improvement of the axial vector currents 

We consider the on-shell 0(a) improvement of the axial vector current both for the heavy- 
heavy and heavy-hght cases. The discussion for the vector case is in parallel with the axial 
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Fig. 3. (a) c^j^^ I c"^ and (b) Cg ^ / c"^ as a function of m^p^ . Circles denote the plaquette gauge action 
and triangles for the Iwasaki gauge action. 



vector case. (See Ref. 1^) for details.) 

The renormaUzed operator with the 0(a) improvement is given by 

A'f'^\x) = Zf; [g(x)7^75Q(x) -/cX9;{g(x)75g(x)} - (7^c:,^9;{g(x)75Q(x)} 

+9'c'^,{diqix)ha^.l5Qix) - 9'c^Jixh^.lM^^Qix)} + Oig^)] , (4-1) 

where we assume that the Euclidean space-time rotational symmetry is not retained on the 
lattice. The coefficients and c^'"'^''^'' are functions of the quark masses tuq and m^. 
With the aid of equation of motion we are allowed to set c^^ = = 0. In the special case of 
tuq = TUq, = and derived from the charge conjugation symmetry. We 

also note that all the improvement coefficients except vanish in the limit of itlq = rriq = 0. 
We determine Z^** and c^a^~'^'^^ at the one-loop level for both heavy-heavy and heavy-light 
cases. In the following we are restricted to the time component of the axial vector current 
because of the limitation of space. As for the space component refer to Ref. Ej). 

The general form of the off-shell vertex function at the one-loop level on the lattice is 
given by 



9, mp,, mp2) = 7o75i^f + 7075^^1" + ^7075^1" + 



05 
4 
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+(Po - go) hGf + i,i>Gf + i^^Gf + 4i^i>Gf] (4-2) 
+ (Po + go) h^f + 751^^^°' + il^Hf + MHf] + 0(a2), 

where p denotes incoming quark momentum and q for outgoing one. The coefficients F"^, 
are functions of p^, g^, p ■ g, m^i and mp2. 
Sandwiching eq. ()4-2|) by the on-shell quark states u{p) and u(g), which satisfy '^u{p) = 
impiu{p) and u{q)4 = imp2u{q), the matrix element is reduced to 

M(g)vlo5 (p, q, rupi, mp2)u{p) 

+ {Po - qo)u{q)u{p) {Gf + irUpiGf + inip^Gf - nipimp^Gf] (4-3) 
+ (Po + qo)u{q)l^u{p) {i/f + zmpiiYf + imp^Hf - mp^nip^Hf] + 0{a^), 

where the coefficients are summarized as 

Xo5 = + ^mpiFf + ^mp2F°5 - mp.mp^Fl^ (4-4) 
>05 = Gf + imp^Gf + zmp2Gf - mp^rUp^Gf , (4-5) 
Zo5 = + impiHf + ^mpsi/g' - nipimp^Hf . (4-6) 

Since the above coefficients contain both the lattice artifacts and the physical contri- 
butions which remain in the continuum, we have to isolate the lattice artifacts in order to 
determine the improvement coefficients in eq. ()4-H) . The improvement coefficients are given 

by 

(4-7) 
(4-8) 
(4-9) 

where the continuum contributions are obtained by employing the naive dimensional regu- 
larization (NDR) with the modified minimal subtraction scheme (MS). Note that c^^ = 
for irtpi = mp2 from the charge conjugation symmetry. 

Here we should remark an important point. (Xqs)''^", (Iqs)^^**, (^os)''^" are process- 
dependent: They are allowed to have different finite values and infrared divergences, which 
are regularized by the fictitious gluon mass in our calculation, for decay and scattering pro- 
cesses. However, once we subtract the continuum counter part, the process dependences 
are canceled out and we are left with finite constants, namely the improvement coefficients 
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' (^05)'=°"^* 




= (^05) 
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(105)^°'^* , 




= (^05' 
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Combining ^7^75 and the wave function renormalization factors, we obtain the renormal- 
ization factor of the axial vector currents: 



7latt 



^cont 



with 



ZQMti^fi) = cosh(mJ,°^) +rtsinh(mj,°^) 
Zq%tti^p2 ) = cosh(mJ,?) + n sinh(mg^) 



(4-10) 



(4-11) 
(4-12) 

(4-13) 



where ^q,^ are found in Ref. 

Employing a set of special momentum assignments p = p* = {pq = impi,pi = 0) and 
q = = {% = if^p2i<ii = 0) or g = = (go = —imp2,qi = 0), where subscripts s and 
d represent the scattering and the decay respectively, we extract the relevant coefficients 
-^05; ^05; -^05 fo^ ^0 ffom the off-shell vertex function ()4-2|) : 



^05 



(4-14) 



iTr 

4 



d d 

7570 + ^mpi— yl{,5^(l + 7o)7fc75 + ^mp2^^!)5 (1 + 7o)7fc75 

Oqk 



dpk 



J p=p*,q=q* 



05) 



-Tr 



d 

^^75(1 + 70) + 2tmpi—A'i^{l + 7o)7fc75 



d N 

05 J 



4 



-Tr 



p=p*,q=qt 



^^5^1 + 70)75 



(4-15) 
(4-16) 



where superscripts s and c? in X05, 105; -^05 represent their momentum assignments, and 
and H*^^ are functions of p^, g^ and p ■ g resulting in 



^^05 




Q^05 




dpi 


p=p*,q=q* 


dqi 


p=p*,q=q* 










dpi 


p=p*,q=q* 


dqi 


p=p*,q=q* 










dpi 


p=p*,q=q* 


dqi 


p=p*,q=q* 



0, 



(4-17) 
(4-18) 
(4-19) 



with 7 = 1, 2, 3. 



The set of improvement coefficients are determined from eqs. ()4-7l - l4-9j) . Assuming that 



the mp2a corrections are negligible, we evaluate the improvement coefficients, except c 



(+,-) 
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2 4 6 8 10 

(0) 

m , a 

pi 

Fig. 4. /A^Q^g , c^' for heavy-light case as a function of m!"^^ . Sohd symbols denote the plaquette 
gauge action and open ones for the Iwasaki gauge action. 



as a function of m^;^ with = 0. In eqs. (j4-15j) . (j4-16j) we find that Iqs ~ ^05 are not 
determined if we set = 0- Therefore one should extrapolate data at non-zero m^g^ to 
'w!p2 = 0- We however keep rrLp2 = 0.0001 in our calculation to determine c^'~'* since the 
difference between the value at rn"p2 = 0.0001 and the one extrapolated to m"p2 = is less 
than 1 %. In Fig. |3] we show numerical results of ^^7075 5 '* for the heavy-light case. 
The solid lines denote the interpolation with a rational expression. The errors are within 
symbols. We can find the heavy-heavy case in Fig. 9 of Ref. Ej). 

We also make a brief comment on our recent work of the 0{a) improvement for the 
heavy-light vector and axial vector currents with relativistic heavy and domain-wall light 
quarks.!^ The most important feature in this calculation is that the renormalization and 
the improvement coefficients of the heavy-light vector current agree with those of the axial 
vector current. We have shown that this is indebted to the exact chiral symmetry for the light 
quark irrespective of the heavy quark action. We have also presented how to implement the 
on-shell improvement on the massive domain-wall quark action, which is required to cancel 
out the infrared divergences generating from the one-loop vertex corrections. 
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§5. Numerical studies in quenched QCD 



Now we have achieved the 0{a) improvement both for the quark action and the axial 
vector current, the next step is to check the effectiveness of the improvement. Our numerical 
studies are carried out in quenched QCD focusing on the restoration of the space-time 
symmetry for the heavy-heavy and heavy-light meson systems. We investigate the dispersion 
relation of moving mesons and the difference of the pseudoscalar meson decay constants 
extracted from the temporal and spatial components of axial vector currents. 

We take the clover action with non-perturbative Cg^P^''^ for light quarks. As for the 
heavy quarks, we replace the massless contribution in Cg "* and c^j^^ of the improved heavy 
quark action with that of the non-perturbative one as 



where the superscript PT represents "perturbative" value. With this replacement 0(a) 
errors are completely removed at mg = 0. This is required by a consistency between the 
light and heavy quarks. We generated 200 configurations with the plaquette and the Iwasaki 
gauge actions on a 24'^ x 48 lattice at a"^ ^ 2 GeV. We employ three light quark masses 
corresponding to mps/my ~ 0.56 — 0.77 and four heavy quark masses covering the charm 
quark mass. For comparison we make another simulation with the clover action both for the 
heavy and light quarks. 

To investigate the dispersion relation of heavy-heavy and heavy-light mesons, it is con- 
venient to define the effective speed of light as 



where E{ps) is the pseudoscalar meson energy with the spatial momentum ps. This quantity 
is supposed to be unity in the continuum limit, which means the restoration of relativistic 
dispersion relation. In Figs. Efa), (b) we plot numerical results as a function of the meson 
mass for heavy- heavy (H-H) and heavy- light (H-L) systems. We observe that around J/lf" 
mass in the heavy-heavy system the deviation from c^s = 1 is equal to or larger than about 
10% for the clover heavy quark action, while less than 4% for the 0{a) improved heavy quark 
action. 

We also measure the space-time asymmetry from the difference of the pseudoscalar me- 
son decay constants determined from the temporal and spatial components of axial vector 



CB = {cl^{mQa)-cl^{0)} + cl 
ce = {cl^iruQa) - cl^{0)} + d 
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Fig. 5. Effective speed of light(top) and pseudoscalar meson decay constants as a function of the 
meson masses for the plaquette and the Iwasaki gauge actions. RHQ denotes the 0(a) improved 
heavy quark action and CL represents the clover action for heavy quarks. 



currents: 

(Ol^flPS) E 
"-'{0|/lf|PS> \p,\' 

which should become unity in the continuum limit. Results are plotted in FigslHtc), (d). We 
find that the asymmetry reaches about 10 ~ 20% for the clover action around J/^ mass in 
the heavy-heavy system , while it is less than 7% for the improved action. These observations 
allow us to conclude that the improved action clearly reduces the errors caused by mqa. 

§6. Conclusion and perspective 

We have proposed a relativistic 0(a) improvement to the heavy quarks on the lattice. 
The idea is based on the relativistic on-shell improvement with the finite rriQa corrections. 
We have shown that the cutoff effects can be reduced to ©((otIqcd)^) putting the {rriQa)"- 
corrections on the renormalization factors of the quark mass and wave function Zq. Our 
relativistic approach has the strong point over the nonrelativistic ones: the finer the lattice 
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spacing becomes, the better the approach works. This is a desirable feature because we can 
take the full advantage of configurations with finer lattice spacing generated to control the 
cutoff effects on the light hadron physics. 

In the next step we have determined the 0{a) improvement coefficients, u, r^, Cb and Ce 
in the quark action up to one-loop level for the various improved gauge actions. While u and 
Ts are determined from the quark propagator, we use the on-shell quark-quark scattering 
amplitude for cb and Cg. The rriQa dependences are examined by making the perturbative 
calculations done in a mga dependent way. Employing the conventional perturbative method 
with the fictitious gluon mass as an infrared regulator we have shown that the parameters 
u, Tg, Cb and in the action are determined free from the infrared divergences. This is 
achieved if and only if the tree-level values for u, r^, Cb and are properly adjusted. 

We have also made the 0(a) improvement of the vector and the axial vector currents up 
to one-loop level. Our calculation is carried out both for the heavy-heavy and the heavy- 
light cases with the various gauge actions. It is explicitly shown that the renormalization 
and improvement coefficients of the heavy-light vector current agree with those of the axial 
vector current, once we impose the exact chiral symmetry for the light quark. 

Given the 0(a) improved quark action and axial vector current, we can check their 
effectiveness by investigating the restoration of the space-time symmetry for the heavy-heavy 
and heavy-light meson systems. We have focused on two quantities: the dispersion relation 
and the difference of the pseudoscalar meson decay constants determined from the temporal 
and spatial components of axial vector currents. Our results show clear improvement for 
both quantities. 

Up to now the 0(a) improvement of the massive Wilson quark action works well: We 
do not encounter any theoretical contradiction in the improvement procedure of the quark 
action and the vector and axial vector currents up to one-loop level, and numerical studies 
clearly show the effectiveness of improvement. These encouraging results lead us to the next 
step. Our ongoing projects are the 0(a) improvement of the four-fermi operators up to one- 
loop level and numerical simulations with two- and three-ffavors of dynamical quarks. We 
also plan a detailed scaling study in quenched QCD and a nonperturbative determination of 
the improvement coefficients. 
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